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ABSTRACT 

This is the third paper in a series which combines A^-body simulations and semi- 
analytic modelling to provide a fully spatially resolved simulation of the galaxy forma- 
tion and clustering processes. Here we extract mock redshift surveys from our simula- 
tions: a Cold Dark Matter model with either flo = I (rCDM) or = 0.3 and A = 0.7 
(ACDM). We compare the mock catalogues with the northern region (CfA2N) of the 
Center for Astrophysics (CfA) Redshift Surveys. We study the properties of galaxy 
groups and clusters identified using standard observational techniques and we study the 
relation of these groups to real virialised systems. Most features of CfA2N groups are 
reproduced quite well by both models with no obvious dependence on Qq. Redshift 
space correlations and pairwise velocities are also similar in the two cosmologies. The 
luminosity functions predicted by our galaxy formation models depend sensitively on 
the treatment of star formation and feedback. For the particular choices of Paper I 
they agree poorly with the CfA survey. To isolate the effect of this discrepancy on our 
mock redshift surveys, we modify galaxy luminosities in our simulations to reproduce 
the CfA luminosity function exactly. This adjustment improves agreement with the 
observed abundance of groups, which depends primarily on the galaxy luminosity den- 
sity, but other statistics, connected more closely with the underlying mass distribution, 
remain unaffected. Regardless of the luminosity function adopted, modest differences 
with observation remain. These can be attributed to the presence of the "Great Wall" 
in the CfA2N. It is unclear whether the greater coherence of the real structure is a 
result of cosmic variance, given the relatively small region studied, or reflects a physical 
deficiency of the models. 

Subject headings: galaxies: clusters: general — galaxies: formation — dark matter — 
large-scale structure of the Universe — methods: miscellaneous 



1. INTRODUCTION 

Groups of galaxies are weak enhancements in the projected galaxy distribution. Velocity 
information from redshift surveys confirms, however, that they are indeed true spatial density 



^diaferio,gamk,jgc,swhite(5)mpa-garching. mpg.de 



- 2 - 



enhancements (e.g. peller &: Huchra 1983 ; Tully 1987; Ramella, Geller & Huchra 1989; Ramella, 
Pisani & Geller 1997), although probably still far from equilibrium (e.g. Nolthenius fc White 1987| ; 
Giuricin et al. 198^ ; piaferio et al. 1993| ; Moore, Frenk & White 1993; [Marnon 1991 ). Groups 



probe the intermediate scale between galaxies and clusters (e.g. Mahdavi et al. 1998| ) and can 
provide important constraints on the amount of dark matter, on the formation and evolution of 





Geller & Huchra 1982; 


Mulchaey 


Ponman, Cannon & Navarro 1998; 



Ramella et al. 1999|) . 



Further information about the dynamics of structure formation come from anisotropics of 
the correlation function in redshift space, which are related to moments of the pairwise velocity 
distribution. The first moment {V12) enters the time evolution equation for the spatial correlation 



function, and has so far been estimated directly only for the IRAS 1.2 Jy survey (Fisher et al. 1994b; 
see also Juszkiewicz, Springel fc Durrer 1999| ; Ferrerira et al. 1999| ). Considerably more attention 
has been devoted to the second moment (T12, which measures the random motions of galaxies and 
can be used to constrain Qq via the cosmic virial theorem (Peebles 1976; Davis &: Peebles 1983; 
Fisher et al. 1994b|; pVIarzke et al. 1995| ). Several authors have shown that (T12 is sensitive to 



the presence of clusters within a survey (pVIo et al. 19931; |Zurek et al. 1994 iMarzke et al. 1995| 



Somerville, Davis fc Primack 1997a ) and so requires a large survey volume to obtain a robust and 



unbiased estimate (Mo et al. 1993; Zurek et al. 1994; Marzke et al. 1995; Somerville et al. 1997a; 



Jing et al. 1998 ). Many early estimates of (T12 were substantially below the values predicted by the 
popular CDM model and its variants (see e.g. Marzke et al. 19"95| and references therein). Larger 
values were found by Mo et al. (1993), by Marzke et al. (1995) for the Center for Astrophysics 
(CfA, hereafter) surveys, and by Jing et al. (1998) for the Las Campanas Redshift Survey (LCRS, 
hereafter) . The latter results are for the largest two surveys currently available and suggest that the 



large velocity bias previously advocated to reconcile CDM models with observations (e.g. Carlberg 
1991; Couchman & Carlberg 1992; Gelb & Bertschinger 1994) may be unnecessary. 



Comparing the observed distribution of galaxies with models of structure formation is prob- 
lematic. In particular, if only the dark matter evolution is simulated directly, as is usually the case, 
then some method must be found to specify the positions, velocities and luminosities of galaxies 



before "mock" redshift surveys can be constructed for direct comparison with real data (e.g. Davis 
let al. 19851 ; [White, Tully fc Davis 1988| ; |Cole et al. 1998|) . Kauffmann et al. (1998; Paper I here- 
after) implemented semi-analytic galaxy formation modelling methods on high resolution A^-body 
simulations of volumes ~ W^h~^ Mpc'^. (Here and below h is the Hubble constant Hq in units 
of 100 km s^^ Mpc~^.) This approach enables us to simulate the physical processes relevant to 
galaxy formation, evolution and clustering, and to explore different assumptions about processes, 
like star formation and feedback, which occur on scales below the resolution limit of our A^-body 
simulations. 



Standard semi-analytic models (White & Frenk 1991; 


Lacey & Silk 1991; 


Kauffmann, White 


& Guiderdoni 1993|; Cole et al. 1994 
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use Press-Schechter theory ( Press &: Schechter 1974| ) and its extensions (Bower 1991; Bond et al. 



19911 ; [Kauffmann fc White 199^ ; [Lacey fc Cole 1993| ) to compute the merging history of dark 



matter halos. Within these halos, simple prescriptions, based on observational data and on more 
detailed simulations, regulate gas cooling rates, star formation rates, stellar population evolution 
and feedback from stellar winds and supernovae. Dynamical friction considerations determine the 
merging rate among galaxies within a common halo, and so the competition between the formation 
of spheroids by merging and the formation of disks through cooling of diffuse gas. Such models do 
not, however, provide detailed information about the clustering and the motions of galaxies. 

In our approach (see Paper I), dark matter halos and their evolution through collapse, accretion 
and merging are followed directly in the A^-body simulations. Other processes are then followed 
as in standard semi-analytic models. Diffuse halo gas is assumed to cool and collect in a disk at 
the centre of each halo. Stars form in such disks according to simple, observationally motivated 
laws. The central galaxy is identified with the most bound particle in the halo. When two or more 
halos merge, the properties of the central galaxy of the most massive progenitor are transferred to 
the most bound particle of the remnant. Other galaxies from the progenitors become "satellites". 
These can merge with the central galaxy of the new halo on a time-scale which is related to 



their progenitor mass as expected from numerical experiments (e.g. [Navarro, Frenk fc White 



1995 ). Galaxy properties such as luminosity, colour, stellar and gas mass, star formation rate and 



morphology evolve according to recipes borrowed from earlier semi-analytic work. We are thus able 



to analyse clustering as a function of galaxy properties and of redshift (Paper I; Kauffmann et al. 



1999b| , Paper II hereafter). 



We investigate two variants of a Cold Dark Matter (CDM) universe: the rCDM model, with 
cosmological density parameter Qq = 1, shape parameter T = 0.21 and Hubble constant Hq = 50 
km s~^ Mpc and the ACDM model, with Qq = 0.3, cosmological constant A = 0.7 and Hq = 70 
km s~^ Mpc As discussed in Paper I, in order to obtain an approximate match between models 
and observations, supernova feedback must be very efficient in the rCDM model, expelling large 
amounts of the reheated gas and so suppressing the formation of galaxies in low-mass field halos. 
On the other hand, feedback must be inefficient in the ACDM models and most reheated gas 
must be retained in order to produce a sufficient number of luminous galaxies. Paper I compared 
the properties of galaxies in our two models with a variety of present-day observational data, 
among these the luminosity function in the B and K bands, the colour distribution, the two-point 
correlation function, the pairwise velocity dispersion, and the mass-to-light ratios of clusters. The 
free parameters which control star-formation and feedback in these models were set by fitting 
the I-band Tully-Fisher relation of Giovanelli et al. (1997) and the mean gas content of a spiral 



galaxy at a circular velocity of 220 km s^^. This "normalisation" choice resulted in fiducial models 
where certain other properties, in particular the luminosity functions, are a relatively poor fit to 
observation. Nevertheless, we used these fiducial models in a follow-up study of the evolution of 
clustering to high redshift (Paper II) and, for consistency, we will continue to use them below, even 
though it turns out that setting parameters to optimise the luminosity functions gives better fits 
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to the abundance of groups in magnitude-limited redshift surveys. 

Here we extract wide-angle mock surveys from the simulations in order to perform a more 
detailed analysis of small-scale galaxy clustering. We compare our models with the northern region 
of the Center for Astrophysics Redshift Survey (CfA2N hereafter; peller &: Huchra 1989 ; Huchra 



et al. 19901 ; |de Lapparent, Geller fc Huchra 199ll ; [Huchra, Geller fc Corwin 1995| ). The CfA2N 



covers a volume comparable to the volume of our simulation box. Moreover, it provides the largest 
published catalogue of galaxy groups. We compare the properties of groups in the CfA2N with 
those of groups extracted from a variety of mock catalogues. In addition, we compare the redshift 
space correlation functions and the pairwise velocity dispersions of the real and artificial redshift 
surveys. 

Because our predicted luminosity functions are sensitive to variations of model parameters 
within their plausible range, it is important to analyse small-scale clustering in ways that are 
insensitive to the exact luminosities assigned to galaxies. To gain more insight into this question, 
we compare results for our fiducial models with results from mock surveys of models in which 
galaxy luminosities are adjusted to reproduce the CfA luminosity function exactly. (The luminosity 
ranking of the simulated galaxies is preserved during this adjustment.) Only the abundance of 
groups and the fraction of galaxies in groups are strongly affected by this change. Even when the 
CfA luminosity function is reproduced exactly, some differences in small-scale clustering remain 
between ACDM and rCDM, and between both and the CfA2N. Such differences may thus provide 
additional constraints on models. 

Our work is the first attempt to compare a wide-angle redshift survey with simulations where 
the physics of the formation and evolution of individual galaxies is treated explicitly. Previous work 
has compared the CfA redshift space correlations and group properties with A'^-body simulations 
either by assigning galaxies to dark matter particles according to some high-peak statistical model, 



or by assigning them to halos based on an assumed halo mass-to-light ratio (Nolthenius & White 



1987; Moore et al. 1993; Frederic 1995b| ; [Nolthenius et al. 1997 ; Somerville, Primack & Nolthenius 



1997b). Some recent studies have attempted to simulate galaxy formation in detail but have been 



forced either to treat volumes which are too small to obtain reliable clustering statistics (e.g. 



Weinberg, Katz fc Hernquist 1998 ; Jenkins et al. 1998 ) or to use a resolution which is too poor to 



follow the formation and evolution of individual galaxies (e.g. [Blanton et al. 199S ; Cen &: Ostriker 



1998 ). In our approach, the physical processes important in galaxy formation are treated in a 
simplified way, but are included ab initio and are followed throughout the evolution of our models. 
Quantities such as the galaxy luminosity function are thus predictions of our scheme, rather than 
being imposed as part of the modelling. In Paper I, we show how differing assumptions about star 
formation and feedback affect predictions for such quantities. 

In Sect. we analyse simulated three-dimensional (3D) groups, defined as sets of galaxies 
each occupying a single dark halo, and we discuss biases in the spatial and kinematic distributions 
of the galaxies relative to those of the dark matter. Sect. then describes how we extract mock 
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redshift surveys from our simulations, while in Sect. |^ we make group catalogues from these mock 
surveys and investigate how their properties and their relation to "real" 3D groups are affected 
by the parameters which define them. In Sect. |5| we compare our simulated groups with groups 
extracted in the same way from the CfA2N, and we study how this comparison is affected by 
differences between the simulated and observed luminosity functions. Finally, in Sect. ^ we 
calculate redshift space correlation functions (,{rp, tt) and pairwise velocity dispersions (Ti2{rp) from 
our mock catalogues and compare with the CfA2N. We conclude in Sect. 0. 



2. GROUPS IN REAL SPACE 

As discussed in Paper I, we identify dark matter halos in our simulations using a friends-of- 
friends groupfinder which links particles closer than 20% of the mean interparticle separation. The 
most bound particle in each halo is taken as its centre, and we calculate a radius r2oo within which 
the average mass density is 200 times the critical density. A 3D group of galaxies is then defined 
as a set of three or more galaxies brighter than some chosen absolute magnitude which lie within 
r2oo of a particular halo centre. 

Groups in the real Universe are identified in redshift surveys. Below, we will compile mock 
surveys from our simulations and catalogue groups identified in redshift space. We will then be 
able to address two issues: (1) how well the physical properties of redshift space groups correspond 
to those of 3D groups; (2) how well the simulation groups reproduce the properties of groups in 
the real Universe. 

In this section, we describe the physical properties of the 3D groups. We compare these with 



the corresponding properties of redshift space groups in Sect. 4.2 



2.1. Galaxy Luminosity Functions within Groups 

We follow galaxy formation and evolution only in dark matter halos containing at least ten 
dark matter particles. This resolution limit implies that our galaxy catalogues are complete to a 
blue band absolute magnitude Mb of about —17.5 + 5Log/i. The luminosity functions of galaxies 
brighter than this limit are shown in Fig. || for our two fiducial models. Solid, long-dashed, dot- 
dashed, and dotted lines are Schechter function fits to the luminosity functions of the CfA Redshift 
Survey (Marzke, Huchra & Geller 1994), LCRS ( |Lin et al. 199^ ), Stromlo-APM ( |Loveday etH 



19921) , and ESQ Slice Project (| ucca et al. 1997[) , respectively. Table || lists the parameters of 



Schechter function fits to the fiducial models and compares them to those of the CfA survey. As 
discussed in Paper I, our decision to normalise the models to the I-band Tully-Fisher relation results 
in both models producing too many bright galaxies. The rCDM model also produces a faint end 
slope which is steeper than the data and is a substantially worse fit to the CfA than ACDM. The 
total blue luminosity density {Lb) = (j)*T{a + 2)L% is 6.6 x IO^/iLqMpc"^ for the rCDM model and 
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1.1 X lO^/iLoMpc"^ for the ACDM model. Marzke et al. (1994) find (2.0 ± 0.9) x W^hLQUpc'^ 
for the CfA sample. We will see that this disagreement between the simulated and the observed 
luminosity functions affects the inferred properties of groups. 

Fig. |l| also shows luminosity functions for galaxies within halos of differing mass M200 within 
the virial radius 7'2oo- Less massive halos (M200 ^ 10^^h~^ Mq, thin solid and dotted lines) con- 
tribute mainly to the faint end of the luminosity function. Bright galaxies come from large halos 
(M200 > W^'^h^^ Mq, long-dashed and short-dashed lines). Note that these halos show a slight 
bump in the luminosity function at Mb — 5Log/i ~ —20.5 and —22, similar to the bump in the 
luminosity function of several real clusters (see e.g. p^rentham 1998| ; Koranyi et al. 1998 ; Molinari 



et al. 199^ ) including Coma ( [Biviano et al. 1995| ). In the models, this effect originates from the 



large merger cross-section of the bright and massive central galaxies which tend preferentially to 
accrete galaxies of somewhat lower mass. 



2.2. Cluster Profiles and Biases 

Dynamical studies of groups and clusters have always produced mass-to-light ratios M/L 



which are substantially smaller than that needed globally to close the Universe (e.g. Car lb erg et| 



al. 1996; Ramella et al. 1997| ). However, it is not clear that the galaxy populations in groups are 



representative of the Universe as a whole (see Paper I). In addition, if galaxies are more strongly 
clustered than the dark matter, multiplying the mean luminosity density by the typical M/L for 
groups and clusters can give a highly biased estimate of Oq. Indeed, previous work analysing mock 
surveys drawn from simulations has suggested that high and low density universes produce group 
catalogues which are not only almost indistinguishable in their observational properties, but also 
quite similar to the CfA groups ( |Nolthenius fc White 198^ ; Moore et al. 1993 ; Nolthenius, Klypin 
& Primack 1997). 

To investigate the distribution of galaxies relative to dark matter in rich clusters, Figs. |2| and 
^ compare profiles of density contrast and velocity dispersion for the dark matter and for galaxies 
of differing luminosity and colour in halos with M200 > lO^^/i~^M0. In order to compute average 
profiles for the whole halo sample, we normalise the positions and the velocities of galaxies and dark 
matter particles to the halo virial radius r2oo and to the circular velocity V200 at ^200 ) respectively. 
We then compute the average profile by superposing all the halos in our sample and giving equal 
weight to each galaxy. Note that each halo contains a galaxy at its centre. If we keep these central 
galaxies when computing the average profile, every halo would contribute a galaxy; we would thus 
obtain a galaxy number overdensity profile several times larger than the dark matter profile at 
radii r < 0.1r200) only because of our arbitrary choice of placing a galaxy at the centre of each 
dark matter halo. In real clusters, this might not be the case. Therefore, we exclude these central 
galaxies when computing the profiles. This exclusion has no effect at the radii r ^ 0.1r2oo plotted 
in Figs. I and |. 
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The upper part of Fig. ^ shows mean number overdensity profiles {6{< r)) for the dark matter 
and for galaxy samples split by luminosity and by colour. The lower part gives the bias b, the 
ratio between the galaxy and the dark matter over densities. In most cases the galaxies are less 
concentrated than the dark matter, resulting in b values less than 1 outside the cluster core. There 
is no strong dependence of bias on luminosity but there is a variation with colour. We divide 
galaxies into a red and a blue sample at the median of the colour distribution. Red galaxies are 



more clustered than blue ones, in agreement with the observed morphology-density relation ([Oemler 



19741 ; [Dressier 1980[ [Postman fc Geller 1984t [Whitmore fc Gilmore l"992[ ; see also Dressier et al 



1997 ). Blue galaxies are strongly anti-biased (6 < 1) in the ACDM model, showing relatively little 
concentration towards the centres of these rich clusters. This is a result of our assumption that 
galaxies lose their gaseous halos, and so their reservoir of new gas, when they fall into a cluster. 
They then redden as their interstellar medium is used up and their star formation rate drops. The 
effect is much weaker in rCDM because infall of galaxies continues at a high rate until z = in 
this model. Notice that the number density profiles for the total galaxy population have identical 
shapes in the two models but are offset in normalisation by about a factor of 3. 

It is interesting to compare these results with observed number density profiles derived 



from the CNOC cluster sample ( Carlberg et al. 1997 ) and the ESO Nearby Abell Cluster Survey, 



ENACS ( [Biviano et al. 1997 ; de Theije Sz Katgert 1998[ ). To obtain surface density profiles from 



the models, we assume spherical symmetry and integrate the three-dimensional profiles along a 
line of sight. At radii ^ 2r200i where the number surface density reaches the background value, we 



approximate the three-dimensional galaxy number density profile with a [Navarro, Frenk fc White 



[(1997) profile. In fact, at these radii, galaxies follow the dark matter distribution closely (Figure 



y). The absolute normalisations are difficult to establish from the observational data. We therefore 
focus on profile shapes and on the relative distribution of galaxies with different properties. In fact, 
the observational surveys find different concentrations between red and blue galaxies (CNOC), and 
between emission-line and non-emission-line galaxies (ENACS). 

In Fig. |3|, we compare these surveys with our models. In the CNOC sample ~ 70% of the 
galaxies belong to the red subsample. In the ENACS sample the non-emission-line and emission- 
line subsamples contain ~ 85% and ~ 15% of the galaxies, respectively. We split our simulated 
galaxy samples into red and blue subsamples in such a way as to reproduce these fractions. We 
normalize both the observed and the model profiles at S(r2oo)- For the ENACS sample, we assume 
''200 = 0.92/i~^ Mpc, the mean value of r2oo for the clusters in both our rCDM and ACDM models. 

These surveys give mean profile shapes very similar to those we find here. In addition, the 
differences between the observed subsamples are intermediate in strength between those predicted 
by our two models. Note however that the results of Fig. |^ are only indicative, because the 
relative distribution of galaxy subsamples is sensitive to the selection criteria. For example, when 
the spectra are used to split the ENACS sample into early-type non-emission-line and late-type 
emission-line galaxies, the difference widens ( [de Theije fc Katgert 1998 ). 
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Finally, Fig. ^ shows, for our models, the velocity dispersion (in units of the circular velocity 
^200) for galaxies within spheres of radius r (in units of r2oo)- These dispersions depend very little 
on the luminosity of the galaxies and are very similar to the dark matter dispersions. There is a 
noticeable colour effect, however, which is quite large in the ACDM model. The velocity dispersions 
of the red galaxies trace those of the dark matter, whereas blue galaxies have dispersions a factor 



1.5 — 2.0 larger. A similar effect is seen in real clusters ( Moss &: Dickens 1977 ; Mohr et al. 1996 



Carlberg et al. 1997 ; Biviano et al. 1997| , 1998); it reflects the fact that the few observed blue 



galaxies have just fallen into the cluster and so are more weakly bound than the red galaxies. 

The weak velocity bias we find for most galaxy samples is consistent with the similarity between 
the moments of the pairwise velocity distribution for the galaxies and for dark matter (see Paper 
I)- 



2.3. Dynamical Properties of 3D Groups 

Table § lists the quartiles of the distributions of a variety of 3D group properties. Groups are 
here defined as halos which contain three or more galaxies brighter than Mb = —17.5 + 5Log/i 
within r2oo- 

The harmonic radius Rh, the one-dimensional velocity dispersion, a, and the total luminosity. 
Lb galaxy groups are computed directly using the positions, velocities and luminosities of the 
galaxies. Assuming virial equilibrium, we can combine these to derive a virial mass estimate, 
Mvir, and so an estimated virial mass-to- light ratio, M^it./Lb- (We review our definition of these 
quantities in the Appendix. Note that Lb includes a correction for galaxies fainter than our 
absolute magnitude limit. On average, ~ 50% and ~ 30% of the total halo luminosity comes from 
this correction factor in the rCDM and ACDM model respectively.) In the simulations the virial 
mass estimate can be compared with the true mass of the group M2oo- 

The first thing to notice from Table |2| is that groups are very similar in the two models. 
They have slightly larger sizes, lower luminosities and larger velocity dispersions in ACDM; these 
differences combine to give typical {M^it./Lb) values which are two times larger in the low density 
model. This surprising result can be attributed primarily to the lower luminosity density in the 
ACDM model. If we multiply the median [M^i^/Lb) of the groups in each simulation by the mean 
luminosity density, we can estimate fio by dividing the result by the appropriate critical density. 
These estimates VLq'^ are given in the last line of Table |2| and in both cases are nearly a factor 
of 2 smaller than the actual value of Q.q. Notice that Mvir actually overestimates the true group 
mass systematically. If we use the true median M200/LB of groups to compute Qf^^ rather than 
our median virial estimate, we get 0.40 and 0.13 for rCDM and ACDM respectively. Clearly, the 
standard assumption that groups have the same mass-to-light ratio value as the universe as a whole 
is untrue in either model. The bias in mass-to-light ratio is actually a function of group mass as 
can be seen in Fig. 15 of Paper I. 
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3. MOCK CATALOGUES 

The simulation boxes have volume ~ 6 x 10^ Mpc^ and ~ 2.8 x W^h^^ Mpc^ for the rCDM 
and the ACDM models, respectively. Both simulations are normalised to match the present day 
cluster abundance. As a result both simulation boxes contain several clusters of mass ~ 10^^h~^ Mq. 
The CfA2N has volume ~ 7 x 10^h~'^ Mpc^ within cz = 12000 km s~^ and contains four Abell 
clusters of richness R = 2 (including the Coma cluster) with masses ~ W^^h~^ Mq. The CfA2N 
covers the right ascension range [8'*, XT'*] and the declination range [8.5°, 44.5°], thereby avoiding 
regions of strong obscuration. It is clearly well suited for comparison with the simulations. 



In Sect. 2.1 we showed that the simulated luminosity functions agree rather poorly with 



the observations. The discrepancy with the CfA luminosity function ( pVIarzke et al. 1994 ) is 



particularly large because the CfA survey is unusually dense: the normalisation <j)* is two times 
larger and the characteristic magnitude is ^ 0.5 mag fainter than the average values derived 
from other surveys. Furthermore, the CfA2N is significantly denser than the CfA survey as a whole, 
presumably because the "Great Wall" dominates this region. 

In order to understand how the luminosity functions of our models affect the properties of 
the groups and the redshift space correlation functions we study below, we have constructed mock 
catalogues from our simulations in two different ways: 

1. We keep the luminosities of the galaxies predicted by the semi-analytic recipes; we refer 
to catalogues made using these luminosities as semi-analytic luminosity function (SALF) 
catalogues; 

2. We alter the luminosities of the galaxies in the simulations in order to obtain an exact match 
to the CfA luminosity function. This is done as follows. For each model we select a set of 
luminosities from the CfA luminosity function constrained so that their number equals the 
number of simulation galaxies with original luminosity Mb < —17.5 + 5Log/i and their sum 
equals the simulated volume times the CfA luminosity density. The new luminosities then 
replace the original ones in such a way that the luminosity ranking of the galaxies is unaltered. 
For each galaxy we compute the shift from the original to the new luminosity. Ninety per cent 
of the luminosity shifts lie in the ranges [0.5,1.0] and [—0.9,0.0] magnitudes in the rCDM 
and ACDM models respectively; the median shifts are 0.8 and —0.5 magnitudes respectively. 
We refer to catalogues constructed using these new luminosities as CfA luminosity function 
(CfALF) catalogues. 

The large scale structure of CfA2N is dominated by the Great Wall and the Coma cluster. The 
center of Coma has celestial coordinates aigso = 12'*57'", ^igso = 28.32° (see e.g. purzadyan fc 



Mazure 1998) and its distance from the Milky Way is ~ 70h ^ Mpc. To compile mock catalogues we 



place a hypothetical observer at a distance d E [68, 72]h ^ Mpc from the most massive cluster within 
the simulation box. The location of the observer is chosen to coincide with the position of a galaxy 
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with properties similar to the Milky Way, i.e. a blue band luminosity Lb ^ 1.9 x 10^'^Lq and mass- 
to-light ratio ~ 5Mq/Lq (e.g. pilmore, King &: van der Kruit 1990| ). If these values are valid for 
a fiducial Hubble parameter Hq = 65 km Mpc~^, we have = —19.3 -|- 5Log/i. We thus look 
for a galaxy with magnitude Mb - 5Log/i G [-18.9,-19.7], stellar mass M £ [4,8] x IO^^H'^Mq, 
and star formation rate SFRg [0.1, lO]Af0 yr^^, typical of a normal spiral galaxy (Kennicutt 1998). 



We also require that this galaxy belongs to a dark matter halo similar to the Local Group halo, 
with total mass M £ [0.1, 1.0] x 10^^/i^^Mq and containing no more than 2 galaxies brighter than 
Mb — 5Log/i = —17.5, the magnitude of M33. We find 32 and 15 galaxies satisfying these criteria 
in the rCDM and ACDM model respectively. 

Once we have found the observer's galaxy, we rotate the reference frame so that the massive 
cluster has the coordinates of Coma. We then compile a catalogue of galaxies with the same 
right ascension and declination range as the CfA2N. To each galaxy we assign the radial velocity 
cz = (v — v/ig) -r/r+r, where v is the galaxy peculiar velocity, Vhg is the observer's peculiar velocity, 
r is the relative position between the galaxy and the observer and is in units of km s^^. The blue 
absolute magnitude Mb yields the apparent magnitude niB = Mb + 25 -|- 5Log(r//i^^Mpc). We 
include all galaxies with cz G [500,15000] km and tub € [10, miim]- Redshift and magnitude 
lower cutoffs avoid including faint objects close to the home galaxy. 

We choose the magnitude limit mum of our mock redshift surveys as follows. The CfA2N 
catalogue has a Zwicky magnitude limit (roughly a i?-band magnitude limit) mnm = 15.5. Because 
the semi-analytic and the CfA survey luminosity densities differ by more than a factor of two, we 
fix m-iim in the SALF catalogues by requiring that the number of galaxies in the mock survey be 
~ 6000, the number of galaxies within CfA2N. We obtain miim ~ 15 for the rCDM model and 
i^iim ~ 16 for ACDM. Errors in magnitude estimates of the CfA2N catalogue and the scatter in 
the correlation between blue and Zwicky magnitudes make a shift in rrin^ by ~ 0.5 mag not at all 



unreasonable (see Marzke et al. 1994 and references therein). The CfALF catalogues have the CfA 



luminosity function and luminosity density by definition, so, in this case, we set rrin^ = 15.5. 

Note that our simulation box is 85 (141) h'^ Mpc on a side for the rCDM (ACDM) model. 
Thus, in order to have a survey with a depth of 15,000 km s~^, one periodic replication of the 
simulation box is required. For = 15.5, galaxies more distant than 85 (141) Mpc are 

always brighter than Mb - 5Log/i = -19.15 (-20.25) for rCDM (ACDM); only the brightest 
galaxies enter the mock catalogue from the replicated regions. Since we concentrate in this paper 
on small scale clustering, we do not expect our results to be seriously affected by this limitation. 

Figs. ^, 1^ and |^ show the CfA2N catalogue and two typical SALF catalogues from the rCDM 
and the ACDM simulations respectively. The rCDM mock catalogue does not look very much like 
the real Universe. There is too much structure within 5000 km and the model fails to produce 
coherent sheets or filaments. These features are common to all our rCDM mock catalogues. The 
ACDM catalogues are in better qualitative agreement with the data. There are large voids and 
filaments extending across the full survey volume. Nevertheless, the structures are not as striking 
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or as sharply defined as in the CfA2N. In addition, there is no structure comparable to the "Great 
Wall". We have searched the simulation boxes for coherent sheets of galaxies and have failed to 
find anything as large as the observed wall (see also pchmalzing 1999| ). It is unclear if this is a 
result of the relatively small volume of our simulations (particularly rCDM) or reflects a significant 
problem for the cosmological models we have studied. CfALF mock redshift surveys are closer in 
appearance to the CfA2N because the change in luminosity function forces a distribution of galaxies 
in redshift closer to that observed (Figure ^). Despite this the qualitative differences in large scale 
structure between the three cases remain. 

In the following sections, we analyse the CfA2N and mock catalogues extracted from the 
simulations using identical techniques. We have compiled an ensemble of ten mock catalogues for 
each simulation in order to assess the robustness of our results. Note, however, that since all ten are 
constructed from the same parent simulation, the scatter between statistics estimated from them 
will underestimate the true sampling variance. 



4. GROUPS IN REDSHIFT SPACE 

A major problem for the study of groups in the real Universe is in finding an objective algorithm 
to define them, given that only partial knowledge of the phase space coordinates is available in real 
galaxy catalogues. This algorithm should select groups which correspond as closely as possible to 
real 3D groups, at least in a statistical sense. In the eighties, two main algorithms were proposed: 
the hierarchical method (e.g. Materne 1975 ; [Tilly 19^ , 1987) and the friends-of-friends algorithm 



(iHuchra fc Geller 1982| ; [Nolthenius fc White 1987|) . 



We identify groups in our mock catalogues with the friends-of-friends algorithm described by 
Ramella et al. (1997). The robustness of methods of this type has been studied by Nolthenius &: 
White (1987), Moore et al. (1993), and Frederic (1995a)| . All of these studies used dark matter- 



only A^-body simulations to study the relation between groups selected from mock observational 
catalogues and genuine virialized systems. Their results were very similar to those we find below - 
with careful parameter choices it is possible to arrange a fair correspondence between the two kinds 
of system and to ensure that the statistical properties of groups estimated from the "observational" 
catalogues are quite similar to those of the "real" systems. 



Sect. 4.1 reviews the friends-of-friends algorithm. In Sect. 4.2, we investigate the dependence 
of the statistical properties of groups on the linking parameters. In Sect. |5|we compare our models 
with the CfA2N groups. We find that variations of the galaxy luminosity function do affect the 
absolute abundances of groups, but have very little effect on their median properties. 
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4.1. Method 



The friends-of-friends algorithm is an approximate method for identifying systems which he 
above a chosen number overdensity threshold. We know galaxy positions in redshift space rather 
than in real space. We therefore need to use two distinct linking lengths, Vq and Dq, for the radial 
velocity coordinate and the coordinates projected onto the sky, respectively. When we know the 
galaxy luminosity function </>(M), Dq defines the number overdensity threshold 



Sn 

n 



47r Z)|] 



A/lin 



(j){M)dM 



-1 



1 



(1) 



where M\\^ = miim — 25 — 5Log(c2;j///o) is the faintest observable absolute magnitude at the fiducial 
redshift czf = 1000 km s^^ within a redshift survey with apparent limiting magnitude mw^. Having 
chosen the threshold 5n/n and the linking length, Vq, we link each pair of galaxies which satisfies 



sin ( ^ ) < DoRij 



\czi — czj\ < VoRij 

where czi and czj are the galaxy radial velocities, 9ij is their angular separation and 



1/3 



(2) 
(3) 

(4) 



where Mij = — 25 — 5Log[(c2;j + czj)/2HQ\. Note that the scaling law in eq. (Q) has been 
questioned by many authors. Specifically, replacing the power 1/3 with 1/2 (see the argument in 
[Nolthenius &: White 1987 ; Magtesyan 1988; Gourgoulhon, Chamaraux &: Fouque 1992| ) drastically 
reduces the correlation between redshift and velocity dispersion observed in the Huchra & Geller 
(1982) group catalogue. However, part of this correlation is related to a selection effect rather than 
to the grouping algorithm, because groups with low velocity dispersion usually have few bright 
galaxies and so can only be seen at low redshift. Here, we use eq. (Q) for consistency with the 
Ramella et al. (1997) catalogue. 

One of the goals of identifying groups in the real Universe is to estimate their number density 
as a function of properties such as luminosity or velocity dispersion. To compute the correct 
abundance of groups, we weight each according to its distance ( [Moore et al. 199^ ). We consider 
groups with > 3 members. We can thus identify a group only when its third-ranked galaxy has 
absolute magnitude Mj < — 25 — 5Log{{cz) / Hq) , where (cz) is the mean redshift of the group. 
Mj determines the radius czj of the sphere within which we could have identified this group. This 
group contributes with weight to the total abundance of groups, where 



3 



CZj 



(5) 
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is the proper volume sampled by the group, to first order in Zj, Q is the solid angle of the survey 
and is the (unknown) cosmological density parameter. We consider groups with (cz) < 12000 
km s~^; therefore the first order correction in the volume is 9% at the most, when Qq = 1. 

Note that the simulation galaxy catalogues are complete to Mb = — 17.5 + 5Log/i. This 
magnitude sets a minimum redshift czmin- groups closer than czmin could contain galaxies fainter 
than Mb in a real magnitude limited survey. Therefore, we consider only redshift space groups 
with (cz) > cZmin both in the mock catalogues and in the CfA2N catalogue. 



4.2. Linking Parameters, Interlopers, and Physical Properties 

From the CfA2N catalogue, Ramella et al (1997) compiled a fiducial group catalogue with 
linking parameters 6n/n = 80, and Vq = 350 km s~^. Here, we show that in the simulations these 
linking parameters give groups in redshift space with similar velocity dispersions and luminosities 
to the 3D groups (Sect 2.3) but with substantially larger sizes. These linking parameters will then 
be adopted for the remainder of our analysis. 

The accidental inclusion of interlopers is one of the major problems in identifying groups in 
redshift space. Even with full knowledge of the galaxy locations in the six-dimensional phase space, 
the concept of interloper is ill-defined because groups are not isolated. Here, we define interlopers 
as follows. Each 3D group contains only galaxies lying within the same dark halo. Each member 
of a redshift space group is associated with some dark halo, but such a group may contain galaxies 
belonging to different dark halos. We define the dark halo containing the greatest number of group 
members as the group halo. All members belonging to another halo are then interlopers, and we 
define the interloper fraction /int as the ratio between the number of interlopers and the total 
number of true members plus interlopers. If all members belong to different dark matter halos, the 
group is spurious and we set /int = 1, because the group halo is undetermined. 

Fig. m shows the dependence of /int on the values of the linking parameters Vq and 5n/n for 
two SALE redshift surveys. All the CfALF and SALE mock catalogues we compiled yield similar 
results. Squares show the median of the distribution of /int; error bars show the first and third 
quartiles. The solid square refers to the fiducial catalogue 5n/n = 80, Vq = 350 km s~^. It is 
apparent that /int reaches a minimum when 6n/n ^ 80 — 100 and Vq is small. Note that this 
minimum value is quite large; more than a third of the assigned members of a typical group are 
interlopers. Fig. |^ also shows that a significant fraction of groups are spurious; for the preferred 
values of 5n/n this fraction is 20 to 30%. We find that ~ 40% of the triplets are spurious, whereas 
only ~ 20% of the groups with four or more members are spurious. This result agrees with the 



suggestion of Ramella et al. (1989) that the physical association of triplets should be considered 
uncertain. 



Fig. 10 shows how the weighted quartiles of the harmonic radius and of the velocity 
dispersion a vary with the linking parameters which define our catalogues. By weighted quartiles 
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we mean the 25, 50 and 75% points of the cumulative distribution when each group is assigned 
weight l/^I'j (see equation ^); these should correspond approximately to the quartiles of a volume 
limited group sample. For comparison we also plot the quartiles of the corresponding distributions 
for 3D groups as listed in Table 2. These figures show that Rh depends only weakly on Vq and is 
a decreasing function of 5n/n; for 6n/n ^ 80 it depends little on either linking parameter. On the 
other hand, a is a strong function of Vq for all 6n/n. This systematic behaviour has been noted 
in all previous investigations of grouping algorithms; see Trasarti-Battistoni (1998) for a recent 
discussion. While our fiducial choice, Vq = 350 km s~^, leads to a a distribution which is similar 
to that of the true 3D groups, taking 6n/n = 80 results in Rf^ distributions which are biased high 
by about a factor of two. 

This bias in group size is related to the rather high interloper fraction noted above. If interlop- 
ers are excluded when calculating group properties, then the remaining galaxies are, by definition, 
all members of the same halo and so give a measure of group size which is statistically similar to 
that found using full 3D information. We demonstrate this in Tables ^ and ^ which compare the 
weighted quartiles of the distributions in mock redshift surveys including and excluding interlopers 
with those of the 3D groups repeated from Table |^. Excluding interlopers has little effect on the 
velocity dispersions of the groups but brings their estimated sizes and mass-to-light ratios into 
much better agreement with those of the 3D groups. Unfortunately, of course, it is not possible to 
exclude interlopers from real redshift survey group catalogues in this fashion. 

It is curious that in our models this interloper-induced overestimation of Rh in redshift survey 



groups combines with the overestimation of true group mass by Mvir (see section 'LS) to almost 
exactly cancel out the factor of 2.5 difference between the true mass-to- light ratio of groups and the 
mean mass-to- light ratio of the Universe. As a result, as may be seen in Tables |^ and ^, estimates 
of r^o obtained by multiplying the (incorrect) mean estimated M^ij-fLs by the average luminosity 
density are fortuitously quite close to the true value. It is difficult to judge whether this interloper- 
bias conspiracy will also work in the real Universe. Some of the tests we carry out below suggest 
that there is no reason to expect that these two biases should always be of the same order. 



5. COMPARISON WITH THE CfA2N GROUPS 

In this section, we compare our groups with the CfA2N catalogue. The galaxy luminosity 
function is a fundamental ingredient in constructing mock catalogues and in the friends-of-friends 
algorithm, so we might expect properties of groups in our simulated redshift surveys to depend 
strongly on the adopted luminosity function. We show below that this is only partly the case. We 
investigate how the luminosity function affects group abundances (Sect. ^^) and the dependence 
of group properties on the linking parameters (Sect. |5.1] ). 
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5.1. Linking Parameters 

In this section we study how the median group velocity dispersion a depends on the hnking 
parameters used to define groups, and we compare with the behaviour seen in the CfA2N survey. 
The harmonic radius and the total group luminosity Ltot are also directly measurable, but, 
as we have shown, the estimated value of is biased high by interloper contamination, while 
Ltot is insensitive to the linking parameters. On the other hand, the velocity dispersion a is 
not significantly biased by interloper contamination and is sensitive to the linking parameter Vq. 
Moreover, the median cr is a direct measure of the dynamics of galaxies on small scales, which in 
turn depends on the underlying mass distribution. 



The left panels of Fig. 11 show that the trend of a with Vq is similar in all our simulated 
redshift surveys and appears to be independent both of cosmological model (ACDM or rCDM) and 
of the adopted luminosity function (SALF or CfALF). The velocity dispersions of the simulated 
groups are systematically larger than the observed values by an amount which varies from about 
10% at low Vo to about 50% at large Vq. This difference is seen in all ten simulated surveys of each 
model. 

Let us define the grouped fraction A'gai/A^tot as the number of galaxies in groups divided by the 
total number of galaxies in the catalogue. Nolthenius et al. (1997) suggest using the dependence 
of the pair (A'gai/A'tot, o') on Vq at fixed 6n/n as a diagnostic to discriminate between models. 
Differences in the small-scale dynamics can show up as a shifting of the tracks in the A^gai/-/Vtot ~ 



a plane. The right panels of Fig. 11 show such tracks for the CfA2N survey and for our simulated 
redshift surveys. Clearly, the choice of the luminosity function affects Ng^i/Ntot strongly but in 
opposite directions for our two models; if the SALF is used, the rCDM catalogue is closer than 
ACDM to the CfA2N; with the CfALF, we obtain the opposite result. Our models yield a fraction 
of galaxies in groups systematically lower than the CfA2N sample. When Vq = 350 km s~^, the 
fraction of galaxies in groups can be from ~ 7% (CfALF ACDM) to ~ 30% (SALF ACDM) smaher 
than in the CfA2N sample. Note however that the CfA2N track is less than two standard deviations 
from the CfALF ACDM track. 



5.2. Group Abundances 



Figs. |l^ and [l^ show the abundance of groups as a function of harmonic radius and velocity 
dispersion derived from our various catalogues. The top panels show results for the semi-analytic 
luminosity function (SALF). The discrepancy with the observations is simply a reflection of the 
overestimate (underestimate) of the total luminosity density of the Universe by rCDM (ACDM) 
as compared with CfA2N. When we impose the CfA survey luminosity function (bottom panels), 
the simulation results are in much better agreement with the observations, especially for ACDM. 



The change in luminosity function from SALF to CfALF has a strong effect on the normalisa- 
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tion of the group abundance function, but rather httle effect on its shape. As a result, the median 
properties of groups do not depend very strongly on the luminosity function. Tables ^ and ^ list 
the weighted quartiles of galaxy properties in the SALF and CfALF catalogues and compare them 
with the CfA2N values. Taking the variations between different mock catalogues into account, the 
median size and velocity dispersion of groups in the ACDM simulation agree reasonably well with 
the observations, if the CfALF is adopted. Moreover, we note that, according to our prescription 
for identifying galaxies in the simulations, the velocity of the central galaxy coincides with that of 
the dark matter particle with the greatest absolute value of the gravitational potential energy. This 
particle can move rapidly whereas the central galaxy is plausibly almost at rest with respect to the 
barycentre of its dark halo. When we replace the central galaxy velocity with the mean velocity 
of the dark particles within r200 5 the median velocity dispersion of groups becomes ~ 10% smaller. 
This effect is sufficient to bring the CfALF ACDM model in good agreement with the CfA2N. On 
the other hand, groups in the rCDM simulation still have velocity dispersions slightly higher than 
observed. 

Note that, on average, the difference between the central galaxy velocity and the barycentric 
velocity of its dark halo is ~ 80 km s~^. Replacing the central galaxy velocity by the barycentric 
motion has a measurable effect on group dispersions because most groups are triplets and quadru- 
plets. The same replacement does not affect the galaxy pairwise velocity statistics computed in 
section |^. 

For both rCDM and ACDM, groups are significantly less luminous than in the real Universe. 
This is a consequence of the presence of the "Great Wall" in the CfA2N region: most of the groups 
lie within this structure at ~ 7000— 10000 km s^^ and so have a median redshift ~ 8000 km s^^. In 
the models, however, the median group redshift is ~ 6000 km s^^, leading to a typical luminosity 
a factor ~ 2 below that found for the CfA2N. 

Note also that although for the SALF catalogues the interloper-bias conspiracy results in 
estimates of $7o which are quite close to the true values, the same conspiracy does not hold for the 
CfALF. Indeed, when our two simulations are forced to have the same luminosity function they 
produce group catalogues which are very similar in most of their properties. In particular, they 
have similar estimated mass-to-light ratios and so lead to similar estimates of Qq despite the factor 
of 3 difference between the true density values. 

6. CORRELATION FUNCTIONS 

To probe dynamics of galaxies on nonlinear scales, knowledge of galaxy peculiar velocities is 
necessary. Even for nearby galaxies direct measurements have ~ 20% uncertainties and only a 
relatively small number of measurements are currently available (see e.g. [Strauss &: Willick 1995 
and references therein). Alternatively, we can study peculiar velocities statistically by using galaxy- 
galaxy redshift space correlation functions. Marzke et al. (1995) compute such correlation functions 
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for the CfA redshift surveys. Here, we analyse our mock catalogues following their procedures which 
we summarise in the next section. 



6.1. Method 



In redshift space, the vector Sj = cziVi locates a galaxy with redshift czi <C c and celestial 
coordinates rj = (ai,6i). For small angular separations (< 50° in our analysis), we define the 
components of the relative separation s = Sj — Sj of a pair of galaxies 



TT 



S l 

IT' 



(6) 



where 1 = (sj + Sj)/2. The two-dimensional redshift space correlation function ^(rp,7r) measures 
the excess probability, compared to a Poisson distribution, that a galaxy pair has separation (rp, vr). 
We estimate ^{rp, vr) by weighting each galaxy according to the minimum variance estimator ( pavis| 
fc Huchra 19821) . 



By inverting w{rp), the projection of ^(rp,7r) onto the axis, we can estimate the spatial 



correlation function ^(r). Assuming a power-law, w{rp) 
A = ro^r(l/2)r[(7-l)/2]/r(7/2). 



Arl^ , gives ^(r) = [r/ro] and 



On non-linear scales, we can then model £,{rp,TT) as (e.g. Fisher 1995) 

/oo 
-oo 



(7) 



where the pairwise velocity distribution /(rp,7r, y) is well approximated by the exponential form 



f{rp,TT,y) = C(cri2)exp 



0'12 



TT -y 



(8) 



where = + y'^, and y is the pair separation in real space along the line of sight. We model 
the mean streaming velocity {vi2{r)) with the scale-invariant solution to the truncated BBGKY 
hierarchy derived by Davis & Peebles (1983), 



Fr 



l + (r/ro)2' 

where the constant F = 1 in the similarity solution. 



(9) 



Eq. (^) is an excellent approximation to the real distribution on very non-linear scales as 
already discussed by 3heth (1996) and Diaferio fc Geller (1996) . Eq. (|^) also describes both dark 
matter particles and galaxies in our models quite accurately ( piaferio et al. 1999| ). On mildly 
non-linear scales, infall skews the distributions. This skewness can be elegantly formalized with an 



Eulerian perturbative approach ( Juszkiewicz, Fisher fc Szapudi 1998 ) 
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On the contrary, Eq. @ is not a good approximation to the mean streaming velocity measured 
in the simulations. In Sect. |6.2| we will see that the uncertainty on i2('^)) does not strongly affect 
the estimate of the pairwise velocity dispersion cri2, provided it is not assumed to be zero on all 
scales. 

In order to estimate the pairwise velocity dispersion ai2{rp) at fixed projected separation rp, 
we first determine the spatial correlation function ^(r); we then perform the integral in equation 
@ by assuming au constant. Note, however, that ai2 does depend on y = (r^ — r^)^/^ because 
the cosmic virial theorem predicts ai2 oc r'^~'^. However, since 7 ~ 2, this procedure yields values 
of CJ12 in reasonable agreement with those computed with the full three-dimensional information. 
Finally, we use the measured ^{rp,iT) to impose the normalisation C{ai2) in equation (|^. This 
procedure decreases the number of degrees of freedom by one. 

To estimate errors on ^(rp,7r) we use a bootstrap procedure with 50 resampled samples. Esti- 
mates of ^{vp, vr) at different separations are correlated. Therefore, determining (T12 by minimizing 
the computed directly with eq. (^) is not the correct procedure. Principal component anal- 
ysis (Kendall 1975, Fisher et al. 1994a) ) transforms a set of correlated quantities {^} into a set 



of uncorrelated quantities {x}. Meaningful results come from the standard analysis applied to 
this latter set. Note that the errors of both the original and the transformed data sets are not 
Gaussian distributed. Therefore, the confidence levels of the computed should be considered 
only indicative. 



6.2. Results 



Fig. 14 shows the ^(rp,7r) maps of our SALE and CfALE catalogues. At rp ^ 5h^^ Mpc, 
i^(rp,7r) in the SALE rCDM catalogue is larger than in the CfALE catalogue; the SALE catalogue 
contains more galaxies in groups and clusters. On the other hand the two ^(rp,7r) maps for the 
ACDM model are very similar, perhaps because the differences between the two luminosity functions 
are smaller. 

We can quantify the differing behaviour of our two models by looking at the parameters tq and 
7 of the real space correlation function derived by fitting a power law to the projected correlation 
function w{rp). Table |^ shows that tq decreases by 35% from the SALE to the CfALE catalogue in 
the rCDM model, while in the ACDM model there is no significant change. None of these models 
matches the CfA2N parameters satisfactorily, perhaps because the CfA2N yields a correlation 
length which is significantly larger than in other surveys: CfA2 South and SSRS2 have ro = 4.75 
and 5.08/i~^ Mpc respectively ( Marzke et al. 1995| ), while the LCRS gives 5.06/i~^ Mpc ( [Jing et 



al. 1998 ); these values are close to those of our models. Note that CfA2 South also has 7 = 1.99, 



so that its correlation function is in close agreement with that of our ACDM model. 



We find a similar result by looking at the projection of £,{rp,TT) onto the vr axis averaged over 
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the Tp interval (r^jr^) 



' p ' p 



(e(7r))=^— ^/ arp,7r)drp (10) 



(Figs. 15 and 16). Independently of the luminosity function, both models predict a correlation 
amplitude somewhat smaller than that of the CfA2N. Typical amplitude ratios between the CfA2N 
and the models are in the range ~ 1.0 — 2.3. These differences are similar to those between the 
CfA2N and CfA2 South or SSRS2 ( [Marzke et al. 1995| ). 

Small variations in (^(vr)) translate into large variations in the best fit parameters ai2- In fact, 
Marzke et al. (1995) show that ai2 can vary by a factor as large as three from survey to survey. 
In Fig. |l7|, solid squares show ai2{rp) when the similarity solution for the mean streaming velocity 
is assumed (F = 1 in eq. |^). The bold solid line is the 0"i2 of galaxies computed with the full 
three-dimensional information. The variations in the amplitude and slope of (C(vr)) between SALF 
and CfALF catalogues make 0"i2 vary by ~ 200 km s~^. However, given the sensitivity of (J12, the 
differences between the models are not significant: both models predict a (7i2{rp) profile close to 
that of the CfA2N, independently of the luminosity function adopted. A more robust statistic, e.g. 



the single-particle- weighted statistic ai suggested by Davis, Miller fc White (1997) , might provide 
a more convincing comparison. 

Finally, we note that in our models {vi2{r)) is significantly different from zero at megaparsec 
scales; moreover, the similarity solution (eq. Q) does not describe its behaviour adequately. How- 
ever, the estimated (T12 are only weakly dependent on the form of {vi2{'f)). For comparison, we 
estimate 0"i2 with F = 2.5 and rg = 3h~^ Mpc which appears to be a better fit to {vi2{r)) in 
our simulations. The agreement between the estimated 0"i2(rp) and the three-dimensional profile 
does not improve substantially (open triangles in the top panels of Fig. 0). On the other hand, 
in agreement with the result of Marzke et al. (1995), assuming F = underestimates cri2{rp) 



significantly (open squares in Fig. 17). This raises a serious question about the Fourier method 



suggested by Landy, Szalay Sz Broadhurst (1998) for estimating ai2{fp)- this method necessarily 



assumes {vi2{r)) = 0. 



7. CONCLUSION 

In Paper I, we simulated the formation, evolution and clustering of galaxies by combining 
dissipationless A^-body simulations and semi-analytic models of galaxy formation. We investigated 
a high density universe (rCDM) and a low density universe with a cosmological constant (ACDM). 
Here, we extract wide angle mock redshift surveys from these simulations and compare them with 
the CfA2N redshift survey by compiling catalogues of galaxy groups and by computing the redshift 
space correlation function ^(rp,7r). 

Despite their very different f^o values, both models yield a reasonable match to the data, 
although both appear slightly less clustered than the CfA2N: for example, ~ 30% (rCDM) and 
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^ 23% (ACDM) of galaxies are in groups rather than the 32% observed; moreover, the amplitude of 

^{rp, vr) is systematically a factor ~ 1.5 — 2.0 smaller than observed. However, this weaker clustering 
is similar to that observed in the southern regions of the CfA redshift surveys and in the LCRS, 
suggesting that galaxies within the CfA2N catalogue might be unusually highly clustered. 

Our semi-analytic modelling predicts luminosity functions which fit the CfA results rather 
poorly. To test whether this disagreement induces clustering differences between our mock surveys 
and the CfA2N galaxy catalogue, we impose the CfA survey luminosity function on the models. 
We find that in fact only the group abundances are strongly affected, becoming substantially closer 
to those observed in the modified catalogues. Quantities connected with the underlying mass 
distribution, such as the median velocity dispersion of groups or the redshift space correlation 
function, show less sensitivity to the adopted luminosity function. 

We compare our models with the CfA2N because (1) this survey provides the largest catalogue 
of galaxy groups currently available, and (2) it covers a volume similar to that we have simulated. 
However, a coherent two-dimensional structure, the "Great Wall", dominates the large scale distri- 
bution of galaxies in the CfA2N and is not reproduced in our simulations. This structure plays a 
significant role in determining the typical luminosities of the observed groups and may also affect 
other group properties. It is unclear whether the apparent deficiency of coherent two-dimensional 
structure in our mock catalogues is a significant problem for our cosmological models, or reflects 
cosmic variance and the relatively small size of our simulation boxes. The modest discrepancies 
between our models and the CfA2N could plausibly disappear when other surveys, covering larger 
volumes and based on better photometric data become available, and can be compared with simu- 
lations of similarly large cosmological volumes. 

Given the limitations of the observational and numerical data we consider in this paper, the 
agreement between theory and observation is remarkable and demonstrates that present day group 
dynamics alone give little information about JIq- Our preference for ACDM over rCDM is based 
on its better agreement with the observed luminosity functions and Tully-Fisher relations (Paper I) 
and with the observed evolution of clustering (Paper II) rather than on any difference in pairwise 
velocities or in the mass-to-light values implied for groups. On the other hand, these clustering 
differences reflect observable differences in the evolution of groups and clusters and of their galaxy 
populations which should be detectable in future surveys to higher redshift. We will investigate 
these issues more thoroughly in a forthcoming paper. 
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A. APPENDIX 



For reference, we review here the formulae we use to compute the physical properties of groups. 
A system with N galaxies has three-dimensional (3D) harmonic radius 



1 



(Al) 



i=i j=i+i 



where Vij are the pairwise galaxy separations. The corresponding harmonic radius for a real system 
at redshift {cz) is 

'n-i n 



Z Hr, 



S ,Si tan(%/2) 



(A2) 



where 6ij are the pairwise galaxy angular separations and Hq is the Hubble constant. 
The one-dimensional velocity dispersion from the 3D velocities is 



a 



whereas from the line-of-sight velocities it is 



^E(v.-(v))^ 



1/2 



(A3) 



N 



,1/2 



— ^(czi - {cz)f 



i=l 



(A4) 



The combination of Rh and a provides the virial mass via the virial theorem 



G 



(A5) 



where G is the gravitational constant. 



The total luminosity of a galaxy group is the sum of the contribution from the N galaxies 
observed and the contribution from galaxies which are too faint to be observed: 



N 



Ltot — -^faint + X/ 



(A6) 



i=l 
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L{a.mt = Volume X (Lfaint) = poo ,/T\jr ^ / L(j){L)dL (A7) 

Jilim nL)dL Jo 



We can estimate the contribution of the faint galaxies by assuming the luminosity function to be 
universal: 

where is the luminosity of the faintest observable galaxy within the group. We adopt the usual 
Schechter form of the luminosity function 

4>{L)dL = (P*{L/L*Y exp{-L/L*)d{L/L*) (A8) 

or equivalently 

cl){M)dM = 0.41n iO(/,nO-°-^(*^-*'^*)("+i) exp[-10-°-^(*^~^^*)](iM (A9) 
when absolute magnitudes M are used. 
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Table 1. Luminosity Function Parameters 





a 


M% - 5Log/i 


</>*/gal mag ^/i^Mpc ^ 


rCDM 


-1.67 ±0.01 


-21.05 ±0.03 


(5.8 ± 1.0) X 10-3 


ACDM 


-1.40 ±0.01 


-20.13 ±0.02 


(4.1 ±0.8) X 10^3 


CfA 


-1.00 ±0.20 


-18.80 ±0.30 


(4.0 ± 1.0) X 10^2 



Note. — Parameters and \-a standard deviations of the Schechter function fit to the model 
luminosity functions (Fig. IT]) and the CfA Redshift Survey (Marzke et al. 1994). 



Table 2. 3D Groups 



rCDM 



ACDM 



iVgal/iVtot 

a 

LogMvir 
M^ir/Msoo 

LogLfi 

Log(Mvir/LB) 

Oest 
"0 



0.30 
0.13/0.21/0.29 
149/213/296 
12.69/13.05/13.47 

0.88/1.37/2.04 
10.53/10.66/10.89 
2.09/2.37/2.61 
0.55 



0.25 
0.15/0.27/0.43 
179/246/335 
12.92/13.30/13.69 

0.76/1.33/2.13 
10.48/10.61/10.83 
2.39/2.67/2.93 
0.18 



Note. — Quartiles of the distributions for the galaxy groups within the box. Rf^, a, Mvir, Lb, 



and M^i^/Lb are in units of h ^ Mpc, km s 
is the dark halo mass within r2oo- 



h ^Mq, h '^Lq, and HMq/Lq, respectively. M200 
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Table 3. Redshift Space Groups: rCDM 



RS RSc 3D 



iVgai/iVtot 0.30 0.20 0.30 

Rh 0.22/0.35/0.55 0.15/0.25/0.32 0.13/0.21/0.29 

a 121/216/336 151/227/378 149/213/296 

LogMvir 12.87/13.25/13.79 12.73/13.15/13.80 12.69/13.05/13.47 

LogLB 10.51/10.69/10.93 10.61/10.77/10.99 10.53/10.66/10.89 

Log(Mvir/LB) 2.14/2.49/2.95 1.92/2.36/2.71 2.09/2.37/2.61 

nf^ 0.72 0.54 0.55 



Note. — Weighted quartiles of the distributions for the galaxy group catalogue from the rCDM 
SALF survey shown in Fig. |6|: RS refers to groups identified in redshift space, RSc to groups 
identified in redshift space with interlopers excluded, and 3D to properties derived from the 3D 
galaxy distribution. Quantities are in the same units as in Table 0. 



Table 4. Redshift Space Groups: ACDM 



RS RS, 3D 



iVgai/iVtot 0.23 0.15 0.25 

Rh 0.35/0.59/0.80 0.18/0.34/0.55 0.15/0.27/0.43 

a 135/222/345 164/223/368 179/246/335 

LogMvir 12.97/13.53/13.91 12.88/13.37/13.80 12.92/13.30/13.69 

LogLij 10.44/10.65/10.90 10.54/10.71/11.00 10.48/10.61/10.83 

Log(Mvir/LB) 2.37/2.82/3.22 2.28/2.64/3.00 2.39/2.67/2.93 

0.25 0.17 0.18 



Note. — Same as Table ^ for the ACDM SALF survey shown in Fig. ^. 
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Table 5. Luminosity Function Effect: rCDM 



SALF 



CfALF 



CfA2N 



A^gal/iVtot 

Rh 
a 

LogMvir 

LogLs 

Log(M^ir/LB) 

Oest 
"0 



0.30 
0.22/0.35/0.55 
121/216/336 
12.87/13.25/13.79 
10.51/10.69/10.93 
2.14/2.49/2.95 
0.72 



0.22 
0.21/0.39/0.59 
134/225/336 
12.80/13.39/13.82 
10.34/10.52/10.79 
2.28/2.76/3.14 
0.43 



0.32 
0.23/0.44/0.71 
99/183/299 
12.60/13.19/13.84 
10.52/10.88/11.19 
1.77/2.43/2.84 
0.20 



Note. — Weighted quartiles of the distributions for the redshift space group catalogues from 
two corresponding rCDM surveys differing in the luminosity function adopted. The SALF survey 
is shown in Fig. ^. CfA2N column is for the groups in the CfA2N catalogue. Quantities are in the 
same units as in Table S. 



Table 6. Luminosity Function Effect: ACDM 



SALF 



CfALF 



CfA2N 



iVgal/iVtot 

Rh 

a 

LogMvir 

LogLfi 

Log(Mvir/LB) 



0.23 
0.35/0.59/0.80 
135/222/345 
12.97/13.53/13.91 
10.44/10.65/10.90 
2.37/2.82/3.22 
0.25 



0.30 
0.31/0.46/0.64 
122/207/305 
12.87/13.38/13.77 
10.45/10.65/10.87 
2.27/2.69/3.10 
0.37 



0.32 
0.23/0.44/0.71 
99/183/299 
12.60/13.19/13.84 
10.52/10.88/11.19 
1.77/2.43/2.84 
0.20 



Note. — 



bame as Table | for two ACDM surveys. The SALF survey is shown in Fig. 0. 



- 30 - 



Table 7. Correlation Function ^(r) 





ro//i~^ Mpc 


7 


Xl8 


rCDM SALF 


5.00 ±0.13 


1.69 ±0.03 


0.69 


rCDM CfALF 


3.20 ±0.28 


1.95 ±0.11 


0.51 


ACDM SALF 


4.45 ±0.11 


2.03 ±0.03 


1.24 


ACDM CfALF 


4.33 ±0.22 


2.02 ± 0.07 


1.24 


CfA2N 


5.59 ±0.24 


1.82 ±0.04 


0.63 



Note. — Parameters of the correlation function fits for individual mock catalogues. The SALF 
catalogues are shown in Figs. ^ and |^. The scatter between values from different mock catalogues 
is about 5%. 
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tCDM 



ACDM 



o 



td 

G 
so 




-23 -22 -21 -20 -19 -18 -23 -22 -21 -20 -19 -18 
Mg-SLogh Mg-SLogh 

Fig. 1. — Upper panels: Blue band galaxy luminosity function for our models. Error bars are 
Poisson 1-a standard deviations. Solid, long-dashed, dot-dashed, and dotted lines are the CfA 
(IMarzke et al. 1994| ), LCRS (|Lin et al. I996D , Stromlo-APM ( |Loveday et al. 1992| ), and ESP 



(Zucca et al. 1997) luminosity functions, respectively. Lower panels: Luminosity function of 
galaxies within halos of different mass: long-dashed line, M200 > W^^h~^MQ] short-dashed line, 
W^'-^h-^MQ < M200 < lO^'^h-^MQ] dotted line, W^^H-'^Mq < M200 < W^^h-^MQ] thin solid line, 
-^200 ^ 10^^/i~^Mq. The bold solid line is the total luminosity function. 



-32- 




Fig. 2. — Average number overdensity profiles {6{< r)) and bias b of galaxies relative to dark 
matter in halos with mass M200 > W^'^h^^ Mq. Central galaxies of each halo arc excluded in the 
computation of the profiles. The two left most (right most) columns are for the rCDM (ACDM) 
model. First and third column: solid, dot-dashed, and long-dashed lines are for dark matter 
and galaxies brighter than Mb — 5Log/i = —17.5, and —19.5, respectively. Second and fourth 
column: solid and dot-dashed (short-dashed) lines are for dark matter and galaxies brighter than 
Mb = —17.5 -I- 5Log/i and with colour B — I < 1.75 {B — I > 1.75), respectively. The colour cut is 
at the median of the colour distribution which is the same for both cosmologies. 
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CNOC 



ENACS 




Fig. 3. — Surface number density profiles of galaxies within halos with mass 

compared with observed cluster samples as explained in the text. Solid (dashed) lines show the 
profiles for red (blue) galaxies in the models. Crosses and triangles arc for red and blue galaxies 
in the CNOC sample (left panels) and for non-cmission-line galaxies and emission-line galaxies for 
the ENACS sample (right panels) . The fraction of galaxies in each subsample is the same for both 
models and observations. 
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Fig. 4. — Velocity dispersion profiles in units of the circular velocity V200 and velocity bias of 
galaxies relative to the dark matter. Panels and lines are as in Fig. 
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CfA3N 




1B50 galasles 30.5" iSA 33.5' 1403 galasles 33.5" iSA 44.5" 



Fig. 5.— 



Galaxy distribution in the CfA2N catalogue projected onto three declination intervals. 
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Fig. 6. — Same as Fig. |5| for a SALF catalogue extracted from the rCDM simulation box. 
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ACDM (SAIF) 




3049 galaxies 30.5* 33.5° 1476 galaxies 32.5" £5i 44.5" 



Fig. 7. — Same as Fig. |5| for a SALF catalogue extracted from the ACDM simulation box. 
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AC DM (.CfALF) 




gl49 gaLaxies 30.5* 33. S"* 1278 galaxies 32.5" £5i 44.5" 



Fig. 8. — The CfALF catalogue corresponding to the SALF catalogue shown in Fig. extracted 
from the ACDM simulation box. 
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Fig. 9. — Upper panels: Median of the fraction of interlopers /int within an individual group in 
group catalogues extracted from the two SALF catalogues shown in Figs. |6| and with different 
linking parameters 6n/n and Vq. Squares are the weighted medians of the distributions, error bars 
show the lower and upper quartiles. At fixed 6n/n, quantities for Vq = 150, 350, and 550 km s~^ 
are shown; quantities for Vq > 150 km are shifted to the right for clarity. Solid square refers to 
the fiducial catalogue. Lower panels: Fraction fsp of spurious groups. 
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Fig. 10. — Dependence of the group median properties on the Unking parameters 6n/n and Vq in 
the SALF catalogues of Figs. |^ and ^. Values of Vq are as in Fig. |9|. Error bars show the upper 
and lower quartiles. Solid lines show the quartiles of the 3D group catalogue. 
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Fig. 11. — Left Panels: Weighted median velocity dispersion a of group catalogues versus the link- 
ing parameter Vq at fixed 6n/n = 80. Triangles, open and solid squares indicate the CfA2N, rCDM 
and ACDM catalogues respectively. Model groups are extracted from the galaxy catalogues with 
the semi-analytic model luminosity function (top panel) or the imposed CfA luminosity function 
(bottom panel). Points show the mean values averaged over the ensamble of ten mock catalogues. 
Error bars are the l-a standard deviations. Solid squares (ACDM model) are slightly shifted to the 
left for clarity. Right Panels: a versus the fraction of galaxies in groups A^gai/A'tot for the same 
catalogues of the left panels. Points from left to right correspond to increasing Vq. Symbols are as 
for the left panels. 
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Fig. 12. — Group abundance by harmonic radius Rh. Number densities are estimated for catalogues 
extracted with 6n/n = 80 and Vq = 350 km s~^. Bold lines are the mean number densities averaged 
over the ensamble of ten mock catalogues. Shaded areas show the S-cr deviations. Dashed lines are 
for the CfA2N groups. Error bars on the CfA2N curves are Poisson 3-a deviations. 
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Fig. 14. — Maps of the two dimensional correlation function ■^(rp,7r). Top panels are for the SALF 
catalogues of Figs. |^ and 0. Bottom panels are for the corresponding CfALF catalogues. The 
bold contour indicates ^,{rp,'K) = 1. Contour levels are separated by logarithmic intervals of 0.1 for 
^(rpjVr) > 1 and by linear intervals of 0.1 for C(rp,7r) < 1. Dashed contours indicate £,{rp,T^) < 0. 
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Fig. 15. — Correlation functions projected onto the it axis for different intervals of the projected 
separation Vp for the SALF rCDM catalogue shown in Fig. ^ (squares). Intervals of Vp in units of 
Mpc are shown in the right upper corner of each panel. Bold solid lines show the best fits to 
eq. ([lO|). Thin solid lines are the best fits for the corresponding CfALF catalogue. Long dashed 
lines are the best fits for the CfA2N catalogue. All curves are for F = 1 in eq. @. {$,{tt)) curves 
for the ten mock catalogues can differ by ~ 10 — 20% depending on the interval. 
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Fig. 17. — Pairwise velocity dispersion ai2 at different projected separations corresponding to 
the curves in Figs. ^ and 16. Solid (open) squares are for F = 1 (F = 0) in eq. @. In the top 
panels, triangles are for F = 2.5 and vq = 3h~^ Mpc in eq. (^. The 0"i2(?'p) profiles have a scatter 
of ~ 15% for different mock catalogues. Bold lines show the three-dimensional pairwise velocity 
dispersion of galaxies brighter than Mb = —17.5 + 5Log/i (Paper I); the dashed line is the CfA2N 
profile when F = I ( Marzke et al. 1995| ). 



